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Abstract. Let {a,k} < j^ 1 be free identically distributed positive non-commuting 
■ random variables with probability measure distribution /i. In this paper we proved 

Ci \ a multiplicative version of the Free Central Limit Theorem. More precisely, let b n = 

& ' aj^a^/ 2 . . . a n . . . a^ 2 a^ 2 then b n is a positive operator with the same moments 



> 



o 



as x n = a\ai...a n and bl/ 2n converges in distribution to positive operator A. 



We completely determined the probability measure distribution v of A from the 
distribution /i. This gives us a natural map Q : A4 + — > with fi i— > 
v. We study how this map behaves with respect to additive and multiplicative 
free convolution. As an interesting consequence of our results, we illustrate the 
| relation between the probability distribution v and the distribution of the Lyapunov 

exponents for the sequence {ak\kLi introduced in [13] . 

+-> 



1. Introduction 

Denote by M. the family of all compactly supported probability measures defined in 
the real line M. We denote by A4+ the set of all measures in M. which are supported 
on [0, oo). On the set M. there are defined two associative composition laws denoted 
by * and EEL The measure is the classical convolution of \x and v. In probabilistic 
terms, fj, * v is the probability distribution of X + Y , where X and Y are commuting 
independent random variables with distributions \x and u, respectively. The measure 
/iEBz/ is the free additive convolution of ji and v introduced by Voiculescu [19] . Thus, 
\i EB v is the probability distribution of X + Y, where X and Y are free random 
variables with distribution \x and is, respectively. 

There is a free analogue of multiplicative convolution also. More precisely, if \l 
and v are measures in Ai + we can define // M v the multiplicative free convolution by 
the probability distribution of X l l 2 Y X 1 ! 2 , where X and Y are free random variables 
with distribution \x and is, respectively. 

In this paper we proved a multiplicative version of the Free Central Limit The- 
orem. To be more precise, let {ak}^ =1 be free positive identically distributed non- 
commutative random variables with distribution fj,. Let us define 
and b n := a\ a^ 2 . . . a n . . . a 2 a*/ 2 . Then x n and b n have the same moments and 
there exists a positive operator A such that 

b l J n — > A 2 in distribution. 
Moreover, if is is the distribution of A, then 

is = f35 + (T with da = f(t)l, 1/2 ± 1/2u ^(t)dt (1.1) 
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where (3 = M{0}), /(*) = (F^ l> )'(t) and = S ^ ~ ^ (^ <_1> is the 



inverse with respect to composition of F^). 

This gives us, naturally, a map 

Q : M+ —> M + with fjLt-^ Q(fi) = u. 

The measure (//) is a compactly supported positive measure with at most one atom 
at zero and {?(/i)({0}) = /i({0}). 

We would like to mention that Vladislav Kargin in Theorem 1 of [12] proved an 
estimate in the norm of the positive operators b n . More precisely, he proved that if 
r(ai) = 1 there exists a positive constant K > such that 



where <r 2 (ai) = r(a 2 ) — r(ai) 2 . 

It is interesting to compare this result with the analogous result in the classical case. 
Let {cik}'j? = i be independent positive identically distributed commutative random 
variables with distribution /j. Applying the Law of the Large Numbers to the random 
variables log(afc), in case log(afc) is integrable, or applying Theorem 5.4 in [6] in the 
general case, we obtain that 



where the convergence is pointwise. 

The Lyapunov exponents of a sequence of random matrices was investigated in the 
pioneering paper of Furstenberg and Kesten [8] and by Oseledec in [17) . Ruelle [T8] 
developed the theory of Lyapunov exponents for random compact linear operators 
acting on a Hilbert space. Newman in jll] and [15] and later Isopi and Newman in 
[TT] studied Lyapunov exponents for random N x N matrices as — > oo. Later on, 
Vladislav Kargin [13] investigated how the concept of Lyapunov exponents can be 
extended to free linear operators (see [13] for a more detailed exposition). 

In our case, given {ak}^ =1 be free positive identically distributed non-commutative 
random variables. Let \i be the spectral probability distribution of a 2 , and assume 
that /i({0}) = 0. Then 



where A is a positive operator. The probability distribution of the Lyapunov expo- 
nents associated to the sequence {a/ c }^ =1 , is the spectral probability distribution 7 
of the selfadjoint operator L := ln(A). Moreover, 7 is absolutely continuous with 
respect to Lebesgue measure and has Radon-Nikodym derivative given by 



\/na(ai) < \\b n \\ < Kn ||ai 





l/2n 



A 



d7(f) = e7(e*)l( 



a 
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(t)dt 
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where the function f(t) is as in equation fll.il) . 

Now we will describe the content of this paper. In section §2, we recall some 
preliminaries as well as some known results and fix the notation. In section §3, we 
prove our main Theorem and study how the map Q behaves with respect to additive 
and multiplicative free convolution. In section §4, we present some examples. Finally, 
in section §5, we derive the probability distribution of the Lyapunov exponents of 
the sequence {a^\^ =x . 



2. Preliminaries and Notation 

We begin with an analytic method for the calculation of multiplicative free con- 
volution discovered by Voiculescu. Denote C the complex plane and set C + = {z e 
C : Im(2:) > 0}, = — C + . For a measure v G M.+ \ {5o} one defines the analytic 
function if) v by 

f 00 zt 

for z £ C\[0, oo). The measure v is completely determined by i\) v . The function 
ip u is univalent in the half-plane iC + , and ^(zC + ) is a region contained in the circle 
with center at —1/2 and radius 1/2. Moreover, ip u (iC + ) R (— oo, 0] = (/? — 1, 0), where 
(3 = ^({0}). If we set Q u = ip u (iC + ), the function ip v has an inverse with respect to 
composition 

Xu ■ Sl v — > iC + . 
Finally, define the S'-transform of v to be 

1 + z 

S v (z) = Xu{z) , z eQ u . 

See [2] for a more detailed exposition. The following is a classical Theorem originally 
proved by Voiculescu and generalized by Bercovici and Voiculescu in [4] for measures 
with unbounded support. 

Theorem 2.1. Let E M + . Then 

S^v{z) = S^S^z) 
for every z in the connected component of the common domain of and S v . 

It was shown by Hari Bercovici in |3] that the additive free convolution of proba- 
bility measures on the real line tend to have a lot fewer atoms. To be more precise. 

Theorem 2.2. Let fi and v be two probability measures supported in R. The number 
a is an atom for the free additive convolution of ji and v if and only if a can be written 
as a = b+c where /i({6}) + z/({c}) > 1. In this case, /lEBz/ ({a}) = /i({6}) + ^({c}) — 1. 
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For measures supported on the positive half-line, an analogous result holds, with 
a difference when zero is an atom. The following Theorem was proved by Serban 
Belinschi in [TJ. 

Theorem 2.3. Let \i and v be two probability measures supported in [0, oo)- 

(1) The following are equivalent 

(a) \i IE v has an atom at a > 

(b) there exists u and v so that uv = a and fi({u}) + f({v}) > 1. 
Moreover, fi({u}) + v({v }) — 1 = IEI v ({a}) . 

(2) ^z/({0}) = max{M{0}),K{0})}. 

In [16] Nica and Speicher introduced the class of i?-diagonal operators in a non- 
commutative C*-probability space. An operator T is i?-diagonal if T has the same 
*-distribution as a product U H where U and H are *-free, U is a Haar unitary, and 
H is positive. 

The next Theorem and Corollary were proved by Uffe Haagerup and Flemming 
Larsen in [9] where they completely characterized the Brown measure of an R— 
diagonal element. 

Theorem 2.4. Let (M, r) be a non- commutative finite von Neumann algebra with 
a faithful trace r. Let u and h be *-free random variables in M, u a Haar unitary, 
h > and assume that the distribution \i h for h is not a Dirac measure. Denote [It 
the Brown measure for T = uh. Then 

(1) is rotation invariant and 

supp(/i T ) = [H/T^li" 1 , \\h\\ 2 ] x p [0,2tt). 

(2) The S -transform of h 2 has an analytic continuation to neighborhood of 
M{0}) - 1,0], Sfc 3 ((/i fc ({0}) - 1,0]) = [IHIzMlfc -1 !!!) and S' h2 < on 

(M{o})-i,o). 

(3) MM) = MM) and M£(0, S h *{t - 1)- 1/2 ) =t forte (^({0}), 1]. 

(4) fix is the only rotation symmetric probability measure satisfying (3). 

Corollary 2.5. With the notation as in the last Theorem we have 

(1) the function F(t) = S h a(t - l)- 1 / 2 : (MM)* 1 ] (ll^llj Ml^lb] ^ 
an analytic continuation to a neighboorhood of its domain and F > on 

(MM),i)- 

(2) [I? has a radial density function f on (0, oo) defined by 

9(s) = ^( F<_1> )'( s ) M(^({o})),F(iM). 
Therefore, \it = /M{0})5 + °" with da = g'(|A|)dm 2 (A). 
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3. Main Results 

Let {ak}^ =1 be a sequence of free positive equally distributed non-commuting op- 
erators. Let us define the operators x n and b n by 

A U V 2 1/2 1/2 1/2 /q \ 

x n := a\a 2 . . . a n and b n :=a{ a 2 ...a n ...a 2 &{ ■ (3.1) 

It is easy to see that b n is positive for all n > 1. In the next Lemma we will show 
that b n and x n have the same moments. 

Lemma 3.1. For all n > 1 and for all k we have that r{x k l ) = t(6^). 

Proof. The case n — 1 it is clear. For n = 2 we have that x 2 = a\a 2 and 6 2 = 
al/ 2 a 2 al/ 2 then 

r (&*) = r^a^a^a^f^a] 12 ) = T(a 1 a 2 (a 1 a 2 ) k ~ 1 ) = t(x 2 ). 
Now we proceed by induction. Assume that T(x k ) = r(b k ) for all k > 1. Let 

j , 1/2 1/2 

a^n+i = a x y n and o n+ i = a{ c n a{ 

where y n = a 2 . . . a n+ i and c n = a} 2 2 a}J 2 . . . a n+1 . . . a^ 2 a 2 . We know by the induc- 
tion hypothesis that r(y™) = t(c™) for all m > 1. Now we observe that 

r(x k n+1 ) = r((a iyn ) k ) (3.2) 

and 

r(b k n+1 ) = r{{a^ 2 c n a^ 2 ) k ) = r^c^c^a 1 ' 2 ) = r((a lC „) fc ). (3.3) 

Finally, using the fact that a\ is free from c n and y n and equations (13. 2p and (13.31) 
we conclude the proof. 

□ 

Now we are ready to prove our main result. 

Theorem 3.2. Let {ak}k be a sequence of free, positive, equally distributed non- 
commuting operators with distribution /i in A4 + . Let b n be as in 113.1}) . The sequence 
of positive operators b\[ 2n converges in distribution to a positive operator A with 
distribution v in A4 + . Moreover, 

u = f3S Q + a with da = f(t) 1 / 1/2 x 1/a i (t) dt 

\\\ a i h 'W a i IbJ 

where (3 = //({0}), f(t) = {F<- l >)'(t) and F,{t) = S,(t - l)' 1 / 2 . 

1 /2 

Proof. Let u a Haar unitary *— free with respect to the family {at]k and let h — % . 
Let T be the i?-diagonal operator defined by T = uh. It is easy to see, by the free- 
ness assumptions, that (T*) n T n and b n have the same distribution. Moreover, by [TU] 

the sequence [(7^*) ri T' ri ] 1//2ri converges in the strong operator topology to a positive 
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operator A. Let v be the probability measure distribution of A. 

If the distribution of is a Dirac delta, p = 8\, then h = vA and 

[{T*) n T n ] 1/2n = [X n (u*) n u n ] 1/2n = y/X. 
Therefore, bl/ 2n has the Dirac delta distribution distribution 5^ and v = 5^. 

If the distribution of is not a Dirac delta, let pt the Brown measure of the operator 
T. By Theorem 2.5 in [10] we know that 

/ \X\ p dp T (X) = lim||T n ||| = limr ( [(T*) n T n ]&) = r(A p ) = / t p du(t). (3.4) 

JC n n n \ J J Q 

We know by Theorem 12.41 and Corollary 12.51 that 

p T = p5 + p with dp(r, 9) = ^-f{r) 1(f^),f„(i)] (r) drdd (3.5) 

where /(*) = [F<- l >)\t) and F^t) = S,(t - l)' 1 ' 2 . 
Hence, using equation (13.4)) we see that 

oo ,2tt ,F„(1) 1 

r p du(r) = / / — r p f{r)drdO = I r p f(r)dr 

JO Jf„(P) 27r ' JFnip) 

for all p > 1. Using the fact that if two compactly supported probability measures 
in A4 + have the same moments then they are equal, we see that 

v = f35 + a with da = fit) 1^(0)^(1)] (t) dt. 



By Corollary 12.51 we know that 

FJl) = 1 1 «i /2 1 1 2 and lim Fjt) = WaV 1 ^]^ 1 
concluding the proof. 

□ 

Note that the last Theorem gives us a map Q : A4 + — > M. + with p \— > Q(p) = v. The 
measure Q{p) is a compactly supported positive measure with at most one atom at 
zero and Q{p)({0}) = p({0}). 

Since 

Q(ji) = f35 + a with da = f(t) 1^(0)^(1)] (*) dt 

and fit) = (F< -1> )'(t) where F^t) = S^t - 1) "V 2 for t G ((3,1]. The function 
S^t — 1) for t G (/3, 1] is analytic and completely determined by p. If p\,pi G A4 + 
and ^(t — 1) = Sfj, 2 (t — 1) in some open interval (a, b) C (0, 1] implies that /ii = /i2- 
Therefore, the map Q is an injection. 
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Remark 3.3. A measure /x in M. + is said M-infinitely divisible if for each n > 1 
there exists a measure \i n in Ai + such that 

H = \i n \i n . . . IEI // n (n times). 

We would like to observe that the image of the map Q is not contained in the set of 
^-infinitely divisible laws since an ^-infinitely divisible law cannot have an atom at 
zero (see Lemma 6.10 in 

The next Theorem investigates how the map Q behaves with respect to additive 
and multiplicative free convolution. 

Theorem 3.4. Let p be a measure in M. + and n > 1. If Q(p) = {35q + a with 
do- = f(t) 1(^09)^(1)] (t) dt then 

g(fx Sn ) = p n 5 + a n with da n = y/nf(t/y/n) l [v ^ ( g„+ n -i } ,y/nF„(i)] (*) dt 
where (3 n = max{0, n(3 — (n — 1)} and 

Q(jP») = (35 + Pn With d Pn = -t^ftf'") l {FM n ,^(l)«](f) dt. 

Proof. Recall the relation between the and transform (see [H]), 

( \<-!> 

By the fundamental property of the i?-transform we have R^mn(z) = nR^(z). There- 
fore, 



\ <-i> 

znRn(z)J = zS fJ m n (z). 



Hence 



thus 



Then 



2; 

—SJz/n) = zS u m n (z) 
n ^ 



S mn(z) = -S^z/n). (3.6) 



it is a direct computation to see that 

F^it) = nF<-^>{t/^) -n + l. (3.7) 
By iterating Theorem 12.21 we see that /i ffln ({0}) = max{0,n/3 — (n — 1)} = (3 n . 



8 GABRIEL H. TUCCI 

Now using Theorem 13.21 we obtain 

£0u ffln ) = p n 5 + a n with da n = y/nf{t/y/n) 1(^ Fm (&±^I) ,^f m (i)](*) dt - 

Now let us prove the multiplicative free convolution part, let then 
Then F Hn(t) = iC(^) and therefore, 



^ 1> (*) = -F?- l> (* 1/B )- (3-8) 

By Theorem 12.31 we now that /z Sn ({0}) = /i({0}) = [3. Therefore, using Theorem 13.21 
again we obtain 

g(fj^) = /35 Q + P n With d Pn = -t^f{t l ' n ) l(F M (/3)»,F M (l)«](t) dt. 

71/ 

□ 

4. Examples 

In this section we present some examples of the image of the map Q. 

Example 4.1. (Projection) Let p be a projection with r(p) = a. Then the spectral 
probability measure of p is p p = (1 — a)5 + a8\. We would like to compute G(p p )- 
Recall that 



Q ! \ Z+ 1 



Therefore, 



F,(t) = { ) and F<- 1 >(t) = T - ¥ . 

Hence, 

G(fip) = (1 - a)5 + a with da = ^ 1 ((WS] (i) eft. 

Example 4.2. Let h be a quarter-circular distributed positive operator, 

dfi h = ^V4 - t 2 l[ 0l 2](£) dt. 
simple computation shows that 

Sh2{z) = ^Ti 

hence by Theorem \3.S\ we see that 

dQ{pffl) = 2t 1 [o,i] 00 dt - 
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Example 4.3. (Marchenko — Pastur distribution) 

Let c > and let fi c be the Marchenko Pastur or Free Poisson distribution given by 

d/j c = max{l - c, 0}S + ^ — — l(a,b){t) dt 
where a = [y/c — l) 2 and b = (-y/c + l) 2 . 

It can be shown (see for example pjj that 

1 



z + c 



Therefore, 
Hence, 



F/j, c (t) = Vt-l + c and F<- l> (t) = t 2 + 1 - c. 



£(/x c ) = max{l - c, 0}5 + o with da = 2t l (v / max{e _ li0} iV5] (*) 



5. LYAPUNOV EXPONENTS OF FREE OPERATORS 

Let {ak}^Li be free positive identically distributed operators. Let \x be the spectral 
probability measure of a\ and assume that /x({0}) = 0. Using Theorem 13.21 we know 
that the sequence of positive operators 



2 

CL1CL2 ■ ■ ■ O^n • • • Oi20j\ 



l/2n 



converges in distribution to a positive operator A with distribution v in M. + . Since 
/i({0}) = 0, this distribution is absolutely continuous with respect to the Lebesgue 
measure and has Radon-Nikodym derivative 

di/(t) = /(t)l (1|ar i|, 2 -i j||ai | |a] (t)dt 

where f(t) = (F<~ 1 >) \t) and F^t) = S,(t - I)- 1 ' 2 . 

Let L be the selfadjoint, possibly unbounded operator, defined by L := ln(A), 
and let 7 be the spectral probability distribution of L. It is a direct calculation 
to see that 7 is absolutely continuous with respect to Lebesgue measure and has 
Radon-Nikodym derivative 



d 7 W=eV(e')l (ln| | a - 1| |- Mn| | ai| | 2] (t) 



dt. 



The probability distribution 7 of L is what is called the distribution of the Lya- 
punov exponents (see [H], [15] and [18] and [13] for a more detailed exposition on 
Lyapunov exponents in the classical and non-classical case). 
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Theorem 5.1. Let {ak}kLi be free positive identically distributed invertible operators. 
Let n be the spectral probability measure of a\. Let 7 be probability distribution of the 
Lyapunov exponents associated to the sequence. Then 7 is absolutely continuous with 
respect to Lebesgue measure and has Radon-Nikodym derivative 

C /7W=eV(e t )l (ln| | arl| |- Mn || ai || 2] (t)^. 

where f(t) = (F^)' \t) and F„(t) = S,(t - l)' 1 / 2 . 

Remark 5.2. Note that if the operators are not invertibles in the \\ ■ || 2 then the 
self adjoint operator L is unbounded. See in the next example the case A = 1. 

The following is an example done previously in [T3] using different techniques. 

Example 5.3. (Marchenko — Pastur distribution) Let {a^.}^ be free positive iden- 
tically distributed operators such that a\ has the Marchenko-Pastur distribution fi of 
parameter A > 1. Then as we saw in the Example \4-3\ in the last section 

du(t)=2tl {VX=IVX] (t)dt. 
Therefore, we see that the probability measure of the Lyapunov exponents is 7 with 

^W = 2e2,1 (iMA-i),Ju.w]W*' 

If X = 1, this law is the exponential law discovered by C.M.Newman as a scaling limit 
of Lyapunov exponents of large random matrices. (See [H], [15] and [H]j. This law 
is often called the "triangle" law since it implies that the exponentials of Lyapunov 
exponents converge to the law whose density is in the form of a triangle. 
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